Superradiance Transition in Transport Through Nanosystems 
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Using an energy-independent non-Hermitian Hamiltonian approach to open systems, we fully 
describe transport through a sequence of potential barriers as external barriers are varied. Analyzing 
the complex eigenvalues of the non-Hermitian Hamiltonian model, a transition to a superradiant 
regime is shown to occur. Transport properties undergo a strong change at the superradiance 
transition, where the transmission is maximized and a drastic change in the structure of resonances 
is demonstrated. Finally, we analyze the effect of the superradiance transition in the Anderson 
localized regime. 
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INTRODUCTION 



Open quantum systems are at the center of many re- 
search fields in physics today, ranging from quantum 
computing to transport in nanoscale and mesoscopic sys- 
tems. In particular, electronic transport in the quantum 
regime can be considered one of the central subjects in 
modern solid state physics P, Transport properties 
depend strongly on the degree of openness of the sys- 
tem. In important applications, the effect of the opening 
is large, and cannot be treated perturbatively. Thus, a 
consistent way to take the effect of the opening into ac- 
count for arbitrary coupling strength between the system 
and the outside world is highly desirable. The effective 
non-Hermitian Hamiltonian approach to open quantum 
systems has been shown to be a very effective tool in 
addressing this issue [3, 0, H, @, 0, H • 

In a typical situation, we have a discrete quantum 
system coupled to an external environment character- 
ized by a continuum of states. Elimination of the con- 
tinuum leads to an effective non-Hermitian Hamilto- 
nian 0, 0, H, H 0] • Analysis of the complex eigenvalues of 
the effective Hamiltonian reveals a general phenomenon, 
namely the segregation of decay widths (corresponding to 
the imaginary part of the complex eigenvalues) . Specifi- 
cally, in a system weakly coupled to the external world, 
all states tend to be similarly affected by the opening, but 
once the coupling reaches a critical value, a sharp recon- 
stitution of the system occurs: almost the entire decay 
width is shared by a few short-lived states, leaving all 
other (long-lived) states effectively decoupled from the 
external world. The analogy between decay width seg- 
regation and Dicke superradiance Q has been pointed 
out in Refs. [J, although Dicke superradiance is as- 
sociated with many-body systems, while width segrega- 
tion occurs also in the one-body case. We will refer to 
this phenomenon as the "superradiance transition" in the 
following. Recently, great attention has been given to 
translating typical quantum optics effects, such as Dicke 
superradiance, into a solid state context [10|. In particu- 
lar, the superradiance effect has been shown to occur in 
several mesoscopic systems [111 ■ 



The effective non-Hermitian Hamiltonian approach to 
open systems has been used mainly under the assump- 
tions of Random Matrix Theory (RMT) [H, More 
realistic systems have also been studied, such as nu- 
clei 14 1 and billiards [l5| . In the latter example, seg- 
regation of resonance widths has already been demon- 
strated experimentally [l6j |. The effective non-Hermitian 
Hamiltonian technique has also been applied to phe- 
nomenological open tight-binding models in solid state 
physics I n t nese papers, the existence of a super- 

radiance transition in such models was shown, but the 
explicit connections to realistic systems were not consid- 
ered. For instance, one might ask whether in a realistic 
situation the coupling to the external environment can be 
increased up to the point where a superradiance transi- 
tion occurs. Also, the energy dependence of the effective 
Hamiltonian is not easy to treat exactly, so one might ask 
in which realistic applications this energy dependence can 
be neglected. 

In this paper, we consider the problem of transport 
through a sequence of potential barriers, see Fig. [TJ 
which can be considered a paradigmatic model in solid 
state physics. This potential profile appears in real ap- 
plications, such as semiconductor superlattices or one- 
dimensional arrays of quantum dots, and has been widely 
discussed in the literature 

The case of equall y sp aced potential barriers has been 
analyzed previously [18j |. Here a different and more gen- 
eral approach to the problem is considered. First, we 
show that for weak tunneling coupling among the wells, 
an energy-independent effective Hamiltonian approach 
produces excellent agreement with an exact (numerical) 
treatment of the problem. Moreover, it is shown that 
even in this simple system a superradiance transition oc- 
curs as the coupling to the external world is increased 
by decreasing the widths of the external potential barri- 
ers. With the aid of the effective Hamiltonian approach, 
we recover several previous results and shed new light on 
the essential features of this well-studied model, allowing 
for a detailed understanding of the resonance structure. 
We emphasize that the powerful effective Hamiltonian 
formalism is not in any way limited to simple models of 
this type, and can be applied to situations where exact 
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treatment is difficult or impossible. In order to show 
this, we also analyze the case of random spacings among 
the potential barriers, and observe the consequences of 
the superradiance transition in the Anderson localization 
regime. 

After briefly reviewing the effective Hamiltonian for- 
malism in Sec. [ill we build the effective Hamiltonian 
model for a sequence of potential barriers in Sec. IIIII 
In Sec. HVl the critical coupling value at which the su- 
perradiance transition occurs is derived, and in Sec. [V] 
we discuss the consequences of this transition on the res- 
onance structure. In Sec. IVI| we show that the maxi- 
mum transmission is achieved at the superradiance tran- 
sition, and we estimate analytically the exponential gain 
in transmission due to the superradiance effect. Finally, 
in Sec. IVIIi we consider the superradiance transition in 
the Anderson localization regime, as a function of the 
disorder strength. 

The effective Hamiltonian approach shows great 
promise in experimental applications, such as quantum 
dots [2lJ and photonic crystals [22| • We also believe that 
the superradiance transition can play a major role in ex- 
plaining many of the results found in open mesoscopic 
systems 2l| . even if this effect has often been neglected 
in the literature. 



II. EFFECTIVE HAMILTONIAN 

We first sketch the essential features of the effective 
Hamiltonian approach to open quantum systems. Details 
of the derivation can be found in Refs. [3. [4l pTT. l23j| . 

Consider a discrete quantum system described by N 
intrinsic basis states \i) coupled to a continuum of states 
|c, E), where c — 1 ... M is a discrete quantum number 
labeling M channels and E is a continuum quantum num- 
ber representing the energy. Let A?(E) be the transition 
amplitude between the intrinsic states and the contin- 
uum. Then the effective Hamiltonian for the intrinsic 
system, which fully takes into account its opening to the 
outside, can be written as: 



H eff (E) = H + A(E) - -W(E) 



(1) 



with 



Wy(E) = 2* Yl AUE)A C 3 {E)\ (2) 

c(opcn) 

where the sum is limited to the open channels, and 

„ f Af(E')A c AE')* 
A. tJ (E) = £ P.v. / dE> A J* > . (3) 

c J 

Assuming Wij(E) and Aij(E) are smooth functions of 
the energy, their energy dependence can be neglected if 
the region of interest is concentrated in a small energy 
window. With the aid of the effective Hamiltonian, the 



transmission T ab (E) from channel a to channel b can be 
determined: 



where 



T ab {E) = \Z a \E)\\ 



N 1 



(4) 



(5) 



is the transmission amplitude. 

We can also write T ab (jE;) in a different way, diagonal- 
izing the effective non-Hcrmitian Hamiltonian H c g. Its 
eigenfunctions |r) and (f| form a bi-orthogonal complete 
set, 

H eS \r)=S r \r), (r\H eS = (f\S;, (6) 
and its eigenvalues are complex energies, 



(7) 



corresponding to resonances centered at E r with widths 
T r . The decay amplitudes Af are transformed according 
to 



^ = ]T^(z|r), j* = Yjr\j)A 



(8) 



and the transition amplitudes are given by 



Z ab (E) 



1 



N 
r—1 



(9) 



The complex eigenvalues £ of H c g coincide with the 
poles of Z(E). It is clear that the properties of the com- 
plex eigenvalues of the effective Hamiltonian are very im- 
portant for understanding the transport properties of the 
system. 

As the coupling between the intrinsic states and the 
external continuum is increased, a rearrangement of the 
widths T r occurs. This rearrangement is usually referred 
to as the "superradiance" transition. 

In order to understand the origin of this transition, 
we can consider a simplified version of Eq. |T]): H e g = 
Ho — Tj'yW, where 7 is a parameter that controls the 
coupling strength with the external world (which now 
we assume to be of the same order of magnitude for all 
the intrinsic states), and H is assumed to be diagonal 
with eigenvalues Eq. For small 7, the first-order com- 
plex eigenvalues of H e g are Si — E^ — \^Wu. If wc 
consider the opposite limit of large 7, Hq can be viewed 
as a perturbation acting on W . Due to the factorized 
structure evident in Eq. ((2]), W has only M non-zero 
eigenvalues for M < N. Thus, only M states will have 
a decay width in the limit of large coupling, while all 
others will have zero width to first order. Therefore, as 
the coupling increases, all widths initially increase pro- 
portionally to 7, but at large couplings only M of the 
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widths continue to increase, while the remaining N — M 
widths approach zero. This simple example suggests that 
a transition between these two regimes may take place at 
a critical value of 7 . Roughly, the transition occurs when 
7/-D f» 1 0, [l3|, 0, [2f| , where D is the mean level spac- 
ing of Hq. Note that the qualitative criterion 7/D « 1 
for the transition to superradiance is valid in the case 
of uniform density of states and negligible energy shift; 
when the density of states is not uniform, the transition 
to superradiance occurs as a hierarchical process [24| . In 
the case of a non-negligible energy shift, see the analysis 
in Sec. HV1 

From the above discussion it should be clear that the 
superradiance transition emerges in the non-Hermitian 
effective Hamiltonian approach as a general phenomenon, 
depending not on the details of the system, but only on 
the factorized structure of W. 



III. EFFECTIVE HAMILTONIAN FOR A 
SEQUENCE OF POTENTIAL BARRIERS 



when distances A, A cxt , and L in Fig. Q] are measured 
in nm (the typical scale in semiconductor supcrlatticcs), 
all energies are calculated in units of 0.038 eV. In the 
following we set L = 2 and Vo = 1000. 

We will now proceed to build an effective non- 
Hermitian Hamiltonian to describe the quantum trans- 
port through a sequence of potential barriers. 

A sequence of N potential wells can be thought of as 
a closed system coupled to the continuum of scattering 
states through the external barriers. Changing the ex- 
ternal barrier widths or heights will change the coupling 
to the continuum. In the limit of low tunneling coupling 
among the wells, the usual tight binding approximation 
can be used to model the closed system: we define the 
intrinsic basis states \i) as the bound states in each poten- 
tial well, corresponding to a certain energy level E Q . Each 
basis state is coupled to its nearest neighbor by the tun- 
neling coupling ft. For small coupling we have [2H l27l|: 

2(T 2 En 

n = 2a\^(x Q )\ 2 = — ° cxp(-aA), (10) 



Let us consider quantum transport through a sequence 
of N+l potential barriers, see Fig.[TJ of width A, height 
Vo, and inter-barrier separation L. The transport prop- 
erties will be analyzed as we change the external barrier 
width A cxt while keeping all the other barriers fixed. 
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FIG. 1: Sequence of potential barriers of finite height and 
width. 

We computed the transmission through this system in 
a standard way, by matching the wave function and its 
derivative in every region, see Fig. [1] Writing the wave 
function in region I as ipi = Be lkx , with k = yJ~E, and in 
region VII as ipvu = Ae lkx +A'e~ lkx , we obtain the trans- 
mission coefficient T(E) = \B/A\ 2 and the reflection co- 
efficient R(E) = \A'/A\ 2 = 1 - T(E). For comparing our 
numerical simulations with experimental results, we note 
that we work in h 2 /2m e = 1 units throughout. Thus, 



where a — y/Vo — Eq, k = \J Eq, ip is a basis wave func- 
tion localized in a single potential well, and xq is a point 
in the middle of a potential barrier immediately adjacent 
to that well. Due to the tunneling coupling among the 
N wells, the eigenenergies of the closed system form a 
miniband around Eq, see Eq. (| L5[) in the next Section. 

The outside world is characterized by the scattering 
states to the left, \L,E), and to the right, \R,E), of the 
sequence of potential barriers. Due to the coupling to the 
scattering states, the states |1) and \N) acquire a finite 
width 7 and an energy shift S, which can be computed 
following Refs. [13, S|- In the case of a varying external 
barrier width A cxt , one obtains: 



8a*E k 



7 - ^(l+aL/2) 



■ exp (-2oA ex t) 



(11) 



5=k 



la k 



-7. 



Note that the shift S vanishes for Eo — Vb/2; otherwise 
the sign of d is given by the sign of Eo — Vq/2. 

Analogous expressions when the external potential 
height I4xt is varied can be computed by extending the 
methods of [27] , but are more complicated and will not 
be reported here. 

We can now write the full effective Hamiltonian for the 
miniband centered at energy Eq as: 



/E Q + 5-^ n 

' n e n 



H, 



off 



n E 
V 00 

Using Eqs. (|3|) and the transmission 



/ 



+0 - 5 7 

(12) 

through the se- 
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quence of potential barriers becomes 



T(E) 



(7/0) 



(13) 



From Eq. (flU)) we see that the spectrum of complex 
eigenvalues £k = E^ — \ Tj. of H c g determines the trans- 
mission through the system. 

In order to show the range of validity of the effec- 
tive Hamiltonian model, we compute the normalized in- 
tegrated transmission: 



S = 



1 



T(E)dE , 



(14) 



where the interval [E m i n , E max ] includes the entire mini- 
band centered at Eq. 

The predictions of Eq. (|13|) are now compared with 
the exact numerical results. The effective Hamiltonian 
approach is expected to break down for small values of 
aA. In Fig. [2j we plot S vs aA for a system of N = 10 
wells, with Eq w 20 fixed, so that a remains constant. 
From the figure we can see that the effective Hamilto- 
nian approximation gives excellent results for a A 1. 
Note also that S is independent of A for A/A oxt = 2 
in the weak coupling limit, a A ^> 1, as indicated by a 
dashed line in Fig. [2] This follows from the fact that 
7/f2 is independent of A when A/A GXt = 2, see Eq. (fTU)) 
and Eq. (TT])- Note that we also compared the results 
obtained with the effective Hamiltonian approach with 
available analytical results found in the literature [l8| . 
and found excellent agreement in the regime a A 1. 



IV. SUPERRADIANCE TRANSITION 

We will now analyze the superradiance transition that 
occurs in the effective Hamiltonian model built in the 
previous section. 

Diagonalization of the intrinsic Hamiltonian leads to 
the energy levels @: 



Wq = Eq ~ 2n cos(irq/{N + 1)) . 



with q = 1 . . . N. Due to coupling with the external 
world, the energy levels w q will acquire decay widths T q . 
These decay widths are the imaginary parts of the eigen- 
values of the effective Hamiltonian, Eq. fT2")) , and, for 
7 <C 1, they can be written as: 



4-, 



N + l 



sm 2 (irq/(N + 1)) 



(16) 



We see that all widths increase proportionally to 7 for 
small coupling. In the opposite limit of large 7, only M 
states (where M is the number of channels) will have a 
width proportional to 7, while the widths of the remain- 
ing states fall off as 1/7, as explained above in Sec. [Til In 
our case we have M — 2, corresponding to one scattering 
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FIG. 2: The integrated transmission S is plotted vs a A for a 
system with N = 10 wells. Different values of aA are obtained 
by varying A while keeping Eq ~ 20 fixed. The symbols in- 
dicate exact numerical results for the sequence of potential 
barriers, while the lines represent the effective Hamiltonian 
prediction. The effective Hamiltonian approximation is ex- 
cellent for a A S> 1. Two values of A/A cxt are shown. The 
quantities plotted are dimensionless. 



channel each on the left and right. The two superradiant 
states correspond to the two non-zero eigenvalues of the 
matrix W (Eq. |1])). 

In order to find the critical value of the parameter 7 
at which the superradiance transition occurs, we may 
analyze the average width (r) of the N — M narrowest 
widths as a function of the coupling 7, Fig. [3l At the 
critical value of 7, the average width (r) peaks and begins 
to decrease. This is the signature of the superradiance 
transition. 

We can evaluate this critical value of 7 using the cri- 
terion discussed earlier in Sec. [H] (T)/D w 1. Con- 
sider first the simpler case of vanishing energy shift S. 
The average width (r) is then given by the pertur- 
bative expression, Eq. (|16p . taking into account that 



(15) (sin 2 (nq/(N + 1))) -> 1/2 for large AT. Moreover, from 



Eq. (fl~5)) we find that for large N the mean level spacing 
becomes D — 4SI/N, so we obtain: 



(T) 
D 



2 n 



(17) 



Thus, for 8 = 0, the criticality criterion (T)/D sw 1 im- 
plies 7 « 2fi. Note that this happens when A cxt = A/2 
(see Eqs. (flOl) and (jTTJ) ) , so that the superradiance tran- 
sition occurs when the external barriers are precisely half 
as wide as the internal ones. 

Typical examples for large and small N are presented 
in Fig. [3l where E = Vq/2 to ensure that the energy 
shift 6 — 0. The N = 100 example in the upper panel 
illustrates that the estimate 7 w 2fl for the critical value 
works very well at large N. 
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FIG. 3: The average width, normalized to the mean level 
spacing, is shown as a function of 7/fi for the case A = 0.15 
and Eq ~ Vo/2. When N is large, the transition to superra- 
diance is shown to occur at 7/Sl sa 2, in agreement with the 
analytical estimation, see text. The solid line corresponds to 
an average over all N widths, while the symbols are obtained 
by averaging over the N — 2 smallest widths. The quantities 
plotted are dimensionless. 

We now turn to the 8^0 case. Both the density of 
states and the resonance widths are modified, as we can 
see using second order perturbation theory in small 7: 

e, = » a+ (.-iV2)^ + (^-7V4-M) £ ^ , 

(18) 

where A qp = (2/7V) 2 (l + (-l) 9+p ) 2 sin 2 4> q sin 2 (f> p and 
4> q = nq/(N+l). 

Clearly, the local level spacings A q (7) = Re(£ q — £ q -i) 
and the local resonance widths ^(7) = — Im(£ q + £ q -i) 
depend on the index q as well as the coupling 7. A rea- 
sonable hypothesis is that the superradiance transition 
occurs when the resonances begin to overlap locally, i.e., 
r q (7) w Dqiff) for some q. We have confirmed numer- 
ically that this local overlap criterion gives an excellent 
approximation for the critical value of 7 at which the su- 
perradiance transition occurs, for any 6. Unfortunately, 
second order perturbation theory does not provide an 
accurate analytical estimate for 7, confirming that the 
physics is highly non-perturbative near the superradiance 
transition. 



V. RESONANCE STRUCTURE 

To show the consequences of the superradiance transi- 
tion on the transport properties, here we analyze the res- 
onance structure, by considering the transmission T(E). 
Note that the resonance structure can be directly re- 
solved experimentally, see [20| . 



In this section, we focus on the case of N = 5 po- 
tential wells. As discussed above in Sec. [TV] a signature 
of the superradiance transition is the segregation of res- 
onance widths above the critical coupling. The system 
under study has two open channels, thus we expect two 
resonance widths (associated with superradiant states) 
to continue increasing above the transition, while the re- 
maining widths approach zero. In Fig. 01 we show the 
trajectories of the complex eigenvalues £ i of the effective 
Hamiltonian as A/A oxt is increased. Note that the real 
parts of the eigenvalues experience a leftward shift with 
increasing coupling, since in this case Eq < Vo/2, so that 
5 < 0. 




2.3175 2.318 2.3185 



E. 
1 

FIG. 4: The evolution in the complex plane of the eigenvalues 
of the effective Hamiltonian is shown as the ratio A/ A cxt is 
varied. A system of 5 intrinsic states in the Eo ~ 2 miniband 
is considered, with A = 0.2. The emergence of two superra- 
diant states is clearly visible above the transition. Here and 
in the following figures, we use units h 2 /2m e = 1, so that 
if lengths A, A ex t, and L axe measured in nm, energies are 
calculated in units of 0.038 eV. 

In Fig. the transmission is shown as a function of 
the energy for several values of A/A cxt . For small cou- 
pling, A/A GXt = 1.5, we have N = 5 narrow resonances 
as expected. As we decrease the external barrier widths, 
the transmission increases and near A/A cxt = 2 two of 
the resonances start to overlap. At A/A cxt = 2.15, they 
merge, forming a broad resonance. This shows that the 
superradiance transition has a clear signature in the reso- 
nance structure. As the external barrier widths continue 
to decrease, the height of the superradiant resonance de- 
creases (see the case A/A ext = 2.4 in Fig. [5]), until it 
disappears entirely for large A/A ext , due to destructive 
interference between the two superradiant states. In this 
limit, N — 2 narrow resonances remains. The fact that 
two resonances disappear for very large coupling is not 
surprising; indeed in the absence of the two external bar- 
riers, we simply have a system of N — 2 wells. What is 
interesting is that the two individual resonances disap- 
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pear long before the external barriers vanish, indeed im- 
mediately after the superradiance transition. Note also 
in Fig.[5]that results obtained from the effective Hamilto- 
nian model are indistinguishable from numerical results 
obtained by matching the wave functions. 

The behavior we have demonstrated for the case of 
N = 5 potential wells generalizes easily to a larger num- 
ber of wells, with important quantitative differences. In- 
deed, for a longer chain, the superradiant states disap- 
pear much faster as we increase A / A ext above the critical 
value, i.e., the superradiance transition becomes increas- 
ingly sharp. Interestingly, the critical value of A/A cxt 
becomes both iV-independent and -Eo-independent in the 
large- N limit, with the transition occurring at A/A ex t = 
2, in agreement with previous results 
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FIG. 5: The transmission is shown as a function of energy 
for a system of TV = 5 potential wells, for the same param- 
eters as in Fig. [4] and several values of A/Aext- The exact 
transmission, black solid line, is compared with the result 
obtained from the effective Hamiltonian model, Eq. (|13[) . in- 
dicated by red circles. The two superradiant states merge at 
A/A GX t = 2.15, and then disappear entirely, leaving behind 
TV — 2 resonances, see A/A cxt = 2.4. 



VI. INTEGRATED TRANSMISSION 



which has the same value of the integrated transmission 
S as the original sequence of TV wells. 



o N=3 
+ N=10 
□ N=50 




A/A 



ext 



FIG. 6: The integrated transmission S is shown as a func- 
tion of A/Aext, for different numbers of wells. The symbols 
refer to an exact numerical calculation, while the effective 
Hamiltonian result is indicated by the black and red dashed 
curves. The horizontal dashed line represents the value of S 
for A = Aext- In this example, we use Eq » 20 and A = 0.15. 
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Another interesting quantity to analyze is the inte- 
grated transmission S, Eq. fT4")) , as a function of A/A ext . 
The quantitative enhancement in S when external bar- 
rier parameters are adjusted is important in applications, 
for instance in the design of electron band-pass filters for 
semiconductor superlattices [18]. From Fig. [SI we see 
that S reaches a maximum as a function of the external 
barrier width. Fig. [5] also shows that the value of S for 
A oxt <§; A is the same as for A cxt = A. This is due to 
the fact that S becomes TV-independent for large TV, as 
shown in Fig. [5] as the external barriers disappear we are 
eventually left with a sequence of TV — 2 potential wells, 



FIG. 7: The critical value of j/fl at which the superradiance 
transition occurs (solid curve) is compared with the critical 
value at which the integrated transmission has a maximum 
(circles). For TV = 100 (upper panel), the two coincide, while 
they differ for TV = 5 (lower panel), see the discussion in the 
text. Here A = 0.2. 

The maximum of the integrated transmission can be 
related to the superradiance transition. In Fig. we 
show the critical value of at which the average width 
(r) has a maximum (signaling the superradiance tran- 
sition, see Fig. [3]), compared with the value of 7/fi at 
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FIG. 8: The transmission gain Smax/Si is shown as a function 
of the internal barrier width A (upper panel) and as a function 
of the internal barrier height Vo (lower panel) . In each panel, 
numerical results obtained for energy Eq ~ 500 (circles) and 
Eo w 2 (squares) are compared with the analytical formula 
(|20l) . indicated by the curves. In the upper panel, we fix 
Vo = 1000, while in the lower panel we fix A = 0.2. All 
data points are for a system of N = 10 wells. The numerical 
prefactor in the analytical result, Eq. (|20[l . is chosen to fit the 
numerical data. 



which the transmission S has a maximum. For large N, 
the transmission maximum is reached precisely at the 
superradiance transition, while for small N the specific 
structure of the resonances influences the exact position 
of the maximum. The relationship between the superra- 
diance transition and the transmission maximum can be 
explained as follows: for small 7, the resonance widths 
increase with 7, and so does the integrated transmission 
S. Once the superradiant states start to form, they in- 
terfere destructively and the associated resonances disap- 
pear, while the widths of the other resonances decrease 
with 7, leading to an overall falloff in S. 

In order to estimate the transmission "gain," we com- 
pute the ratio of the maximum transmission to the trans- 
mission at A/A oxt = 1: S max /Si. The transmission is 
given by the area under the resonances, and for isolated 
resonances we can estimate S « 27V(r)/4fi. Since for 
small 7, see Eq. (fTo)) . we have (T) = 2j/(N + 1), we can 
write: 



7 



N 



riN + i 



2 



(19) 



which is iV-indepcndent for large N as noted above. Tak- 
ing into account that the superradiance transition occurs 
at 7 = 2Q, for Eq = Vq/2, the gain can be estimated as: 



o 



(20) 



Si 71 Aol-JEq~ 

In Fig. [51 we show that our estimate works very well 
Eq. (|2U)1 is in agreement with the results obtained in [l8| 



using a different approach. Note also that the gain is 
exponential in the internal barrier width A. 



VII. ANDERSON LOCALIZATION REGIME 

In the previous sections we considered an effective 
Hamiltonian (fT2"]) with equal diagonal energies Eq. Here 
we want to apply the effective Hamiltonian technique to 
the case of random variations of the diagonal energies: 
Eq ± SEq, where SEq is a random variable uniformly 
distributed in [— W/2, +W/2], and TV is a disorder pa- 
rameter. A first analysis of this model can be found in 
Ref. [H. 

Random variation in the diagonal energies can be 
thought of as a consequence of small random fluctuations 
5L of the well widths L. For Eq <C Vo, the eigenenergies 
of a finite potential well may be approximated by the 
eigenenergies of an infinite potential well, Eq — n 2 ir 2 /L 2 , 
where n = 1,2, .... For small fluctuations SL/L <C 1, we 
have 



2n 2 n 2 r 
SEq = -j^SL 



-C5L, 



(21) 



where C — 2n 2 7r 2 /L 3 . Thus, a random variation of 8Eq 
in [— W/2, +W/2] corresponds to a random variation of 
SL in [-W/2C,+W/2C\. 

The effective non-Hermitian Hamiltonian with diago- 
nal disorder is equivalent to an open Anderson tight bind- 
ing model [1, H(| • The eigenstates of the Anderson model 
are exponentially localized on the system sites, with ex- 
ponential tails given by exp(— x/L\ oc ), where for weak 
disorder, the localization length £j oc at the center of the 
energy band can be written as [311 ]: 



'W 

U, = :i nr,.2 ( — 



(22) 



For L\ oc <C N, the transmission decays exponen- 
tially with N; this is the localized regime. Note that 
for zero disorder, the transmission is N- independent, 
as we showed in the previous section. The condition 
L\ oc = N defines a critical value of (W/ft) CI for the local- 
ized regime, at any given N. In the localized regime, the 
transmission is log-normally distributed, and we have [l[ : 



InT) 



N 

2- h Const. 

E\ oc 



(23) 



In Fig. [HI we show the average transmission as a func- 
tion of the disorder strength for two different degrees of 
opening of the system. The results obtained using an ef- 
fective Hamiltonian with diagonal disorder are compared 
with numerical simulations for the disordered sequence of 
potential wells. The agreement is excellent up to a large 
value of the disorder, where of course our approximations 
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FIG. 9: The average of the logarithm of the transmission is 
plotted as a function of W/Q for A/A ex t = 1 and A/A cx t = 2. 
Results obtained from the effective Hamiltonian (curves) are 
compared with results obtained numerically for the sequence 
of potential barriers (symbols). Here we take iV = 50, Eq ~ 
20, and A = 0.2. The dotted vertical line indicates the critical 
value of W/Q, obtained from the condition L\ oc — N, where 
iioc is given by Eq. (f22l) . 



break down. [Indeed O, Eq. (fTTJ|) , has been computed as- 
suming that the energy levels are aligned, which is not 
true anymore in the presence of disorder.] We stress that 
Eq. Q13[) for the transmission remains valid even in the 
disordered case, making the use of the effective Hamil- 
tonian formalism very efficient, since only the eigenval- 
ues of the effective Hamiltonian are needed to obtain the 
transmission. 

The phenomenon of Anderson localization was studied 
in a closed disordered chain or for fixed opening, while 
in our case we can vary the degree of openness of the 
system. The effect of the opening on Anderson localiza- 
tion is not obvious. Will a maximum of the transmission 
still exist as we vary the coupling of the system with the 
external world? Will the localization length change as 
we open up the system? To answer these questions, we 
have analyzed the effective Hamiltonian, neglecting the 
role of the energy shift, i.e., we set <5 = in the following. 
In Fig. [TH1 we compute the average of — (lnT) over 10 5 
realizations as a function of 7/O. The energy is fixed 
at E — Eq. As we vary 7/O, the average transmission 
reaches a maximum, just as in the disorder-free case. 

Interestingly, as the disorder strength increases, the 
transmission maximum (associated with the superradi- 
ance transition) shifts to ever higher values of the cou- 
pling strength 7. Indeed, the mean level spacing D 
increases with growing disorder, so that the condition 
(r)/D w 1 for the superradiance transition to occur 
will be satisfied at increasingly larger values of 7/O. 
For weak disorder, the disorder-induced correction to D 
is second order in the disorder strength, so (j/Q) cr = 



2 + 0(W 2 /Q 2 ). In the opposite regime of large W/£l we 
find D ps 2Q(1 + W~/2ti)/N, and the critical coupling 
is predicted to be (7/f2) cr sa 1 + W/2Q. This estimate 
works quite well, as shown by the vertical lines in Fig. 1101 
Note that the curves shown in Fig. [10] have been found 
to be independent of N. 
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FIG. 10: The mean logarithm of the transmission is plot- 
ted as a function of 'j/Q, for different values of the disorder 
parameter W/Q. Here we take N = 100, -Bo ~ Vb/2, and 
A = 0.6. The dashed vertical lines indicate the predicted 
value of y/Q at which the transmission maximum is expected 
for W/Q = 2, 10, and 20, see the discussion in the text. Note 
that (lnT)o stands for the value of (lnT) at A = A ex t. 

Finally in Fig.QT] we show (lnT) versus TV for W/Q, = 
2, and for two different values of the external barrier 
width: A/A cxt = 1 and A/A cxt = 2, where the max- 
imum of the transmission occurs. Fig. [11] shows that 
Eq. (|23[) works very well in both situations, even though 
the transmission is enhanced when A/A cxt = 2, for all 
values of N. Thus, the localization length in a disor- 
dered ID model is not affected by the opening, but the 
transmission is. 



VIII. CONCLUSION 

We have analyzed quantum transport through a finite 
sequence of potential barriers, a paradigmatic model for 
transport in solid state physics. In this paper, the ef- 
fective non-Hermitian Hamiltonian approach has been 
used to analyze the transmission through this class of 
systems. The main results of our work arc the following: 
i) we show that for weak or moderate tunneling coupling 
among the potential wells, the system is well described 
by an energy-independent effective Hamiltonian. Knowl- 
edge of the complex eigenvalues of the effective Hamil- 
tonian is sufficient to study transport through the sys- 
tem, ii) As the coupling to the continuum is increased 
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FIG. 11: The mean logarithm of the transmission is plotted as 
a function of TV for A/A cx t = 1 (black circles) and A/A ex t = 
2 (red squares) for the case W/Q. = 2, Eo « Vb /2 and A = 0.6. 
The remaining parameters are the same as in Fig. [5] From 
the fit we have (InT) = a - vN, with v » 0.078, for both 
A/Aext = 1, 2, in good agreement with the theoretical result 
v = 2/L loc = 0.076, see Eq. J23). 



this transition, m) The superradiance transition has 
strong effects on the transport properties: specifically, 
the transmission through the system is maximized at the 
superradiance transition. An expression for the transmis- 
sion gain due to the superradiance transition is derived. 
Moreover the resonance structure is drastically affected: 
at the superradiance transition, we have the formation 
of a broad resonance corresponding to the superradiant 
states. Beyond the transition, this broad resonance dis- 
appears, and the number of resonances decreases by two. 
iv) The case of a disordered sequence of potential barri- 
ers has been also analyzed. In the presence of disorder, 
Anderson localization occurs. We have shown the local- 
ization length remains constant as the opening changes, 
but the transmission has a maximum as a function of 
the coupling to the external world. The critical value of 
the coupling increases with the degree of disorder, and 
we obtain an estimate of the critical value for strong and 
weak disorder, based on the superradiance mechanism. 

In the future it will be interesting to study the conse- 
quences of the superradiance transition beyond the single 
particle approximation, where electron-electron interac- 
tions play an important role. 
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